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Abstract 

We present our study on atomic, electronic, magnetic and phonon properties of one di- 
mensional honeycomb structure of molybdenum disulfide (M0S2) using first-principles plane 
wave method. Calculated phonon frequencies of bare armchair nanoribbon reveal the fourth 
acoustic branch and indicate the stability. Force constant and in-plane stiffness calculated 
in the harmonic elastic deformation range signify that the M0S2 nanoribbons are stiff quasi 
one dimensional structures, but not as strong as graphene and BN nanoribbons. Bare M0S2 
armchair nanoribbons are nonmagnetic, direct band gap semiconductors. Bare zigzag M0S2 
nanoribbons become half-metallic as a result of the (2x1) reconstruction of edge atoms and 
are semiconductor for minority spins, but metallic for the majority spins. Their magnetic mo- 
ments and spin-polarizations at the Fermi level are reduced as a result of the passivation of 
edge atoms by hydrogen. The functionalization of M0S2 nanoribbons by adatom adsorption 
and vacancy defect creation are also studied. The nonmagnetic armchair nanoribbons attain 
net magnetic moment depending on where the foreign atoms are adsorbed and what kind of 
vacancy defect is created. The magnetization of zigzag nanoribbons due to the edge states is 
suppressed in the presence of vacancy defects. 

Introduction 

Unique honeycomb orbital symmetry underlies the unusual properties of two dimensional (2D) 
single hexagonal structures, such as graphene,- - - silicene and group III-IV binary compounds. 4 
Moreover, quasi ID nanoribbons and flakes of these 2D layers have added interesting electronic 
and magnetic properties, which are expected to give rise to important future applications in nan- 
otechnology. 5 ~ 9 Recently, 2D suspended single layer molybdenum disulfide, M0S2 sheets with 
honeycomb structure have been produced.— ^ Single layer M0S2 nanocrystals of ~ 30 A width 
were also synthesized on the Au(l 11) surface and its direct real space STM have been reported. 12 
Unlike graphite and hexagonal BN, the layers of M0S2 are made of hexagons with Mo and S2 
atoms situated at alternating corners. Apparently, 3D graphitic bulk structure called 2H-M0S2, 2D 
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single layer called IH-M0S2, quasi ID nanotubes and nanoribbons of M0S2 share the honey- 
comb structure and are expected to display interesting dimensionality effects. 

Properties of M0S2 nanocrystals are explored in diverse fields, such as nanotribology, 14 ' 15 
hydrogen production, 16,17 hydrodesulfurization catalyst used for removing sulfur compounds from 
oil,— - — solar cells,- 25 photocatalysis. 26 Triangular M0S2 nanocrystals were obtained as a function 
of size by using atom-resolved scanning tunneling microscopy. 27 Photoluminescence emerging 
from IH-M0S2 was observed. 28 Superlow friction between surfaces coated with IH-M0S2 has 
been measured much recently. 15 Using electrochemical methods micro and nanoribbons have been 
synthesized from crystalline 2H-M0S2. 29 Various properties of 2H-M0S2 (see Ref.[— "^ 5 -]), 1H- 
M0S2 (see Ref.[— i 4 - 4 -! 4 - 6 - - — ]) and its nanoribbons (see Ref.[— i 4 - 8 -]) have been an active subject of 
theoretical studies. 

In this paper, we present our systematic theoretical investigation of optimized atomic structure 
and phonon spectrum, mechanical, electronic, magnetic properties of armchair (A-M0S2NR) and 
zigzag (Z-M0S2NR) nanoribbons. Our study reveals interesting results, which are important for 
further study and applications of these nanoribbons. These are: (i) We demonstrated the stability 
of M0S2 nanoribbons through first-principles calculations of phonon frequencies. Specifically we 
deduced the branch of twisting modes, (ii) We calculated force constants and in-plane stiffness 
of armchair and zigzag nanoribbons showing that they are stiff materials, (iii) We examined the 
effects of the reconstruction of edge atoms and their passivation by hydrogen on the electronic and 
magnetic properties of nanoribbons. The energy is optimized through a (2 x 1) reconstruction of 
edge atoms of zigzag nanoribbons, that, in turn, renders half-metallicity. (iv) The properties of 
M0S2 nanoribbons can be dramatically modified by foreign atom adsorption and vacancy defects. 
Since recent works^ 5 -^" 5 -^- show that single layer M0S2 flakes as large as 200 /im 2 can now be 
produced and also can be characterized, present results are crucial for further research on M0S2 
nanoribbons. 
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Method 



Our results are based on first-principles plane wave calculations within density functional theory 
(DFT) using projector augmented wave (PAW) potentials. 52 The exchange correlation potential 
is approximated by generalized gradient approximation (GGA) using PW91 53 functional both for 
spin-polarized and spin-unpolarized cases. All structures are treated using the periodic boundary 
conditions. Kinetic energy cutoff, Brillouin zone (BZ) sampling are determined after extensive 
convergence analysis. A large spacing of ~ 10 A between the S planes of two M0S2 layer are 
taken to prevent interactions. A plane- wave basis set with kinetic energy cutoff of 600 eV is 
used. In the self-consistent field potential and total energy calculations BZ is sampled by special 
k-points by using Monkhorst-Pack scheme. 54 For nanoribbons, BZ is sampled by 1x1x9 k-points. 
All atomic positions and lattice constants are optimized by using the conjugate gradient method, 
where the total energy and atomic forces are minimized. The convergence for energy is chosen as 
10~ 5 eV between two consecutive steps, and the maximum Hellmann-Feynman forces acting on 
each atom is less than 0.05 eV/A upon ionic relaxation. The pressure in the unit cell is kept below 1 
kBar. The phonon dispersion curves are calculated along symmetry directions of BZ within density 
functional theory using Small Displacement Method (SDM). 55 Numerical calculations have been 
performed by using VASP. 56:57 Bader analysis is used for calculating the charge on adatoms. 58 



Properties of Two Dimensional M0S2 

For the sake of comparison we first present a brief discussion of the properties of 2D IH-M0S2 
calculated with the same parameters used for quasi ID nanoribbons. Single layer M0S2 structure 
consists of monatomic Mo plane having a 2D hexagonal lattice, which is sandwiched between two 
monatomic S planes having the same 2D hexagonal lattice. Mo and S2 occupy alternating corners 
of hexagons of honeycomb structure. The contour plots of calculated charge density and difference 
charge density isosurfaces clarify the charge distribution in layers of M0S2 structure. Electronic 
charge transferred from Mo to S atoms gives rise to an excess charge of 0.205 electrons around 
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each S atom. " 6U This situation implies that IH-M0S2 can be viewed a positively charged Mo 
planes between two negatively charged planes of S atoms and this is the main reason why flakes 
of M0S2 structure are good lubricant. The cohesive energy of IH-M0S2 is calculated as 15.55 eV 
using GGA+PAW. The structure is optimized to yield hexagonal lattice constant, a =3.20 A, and 
internal structure parameters, such as the bond distance between Mo and S atoms Jmo-5=2.42 A, 
the distance between two S atoms at each corner ds_s=3.13 A, and the angle between Mo-S bonds 
®5-Mo-5=80.69°. The ground state of monolayer IH-M0S2 is nonmagnetic semiconductor having 
direct band gap of £g=1.58 eV. The upper part of the valence and the lower part of the conduction 
bands are dominated from bonding and antibonding Mo-Ad and S-3p orbitals. 

ID MoS 2 Nanoribbons 

Two dimensional IH-M0S2 can maintain its physical properties, when its size is large. However a 
small flake or a ribbon can display rather different electronic and magnetic properties. In particu- 
lar, edge atoms may influence the physical properties. The passivation of edge atoms by hydrogen 
atoms also result in significant changes in the properties of the nanoribbons. In this respect, one ex- 
pects that the armchair (A-M0S2NR) or zigzag (Z-M0S2NR) nanoribbons of IH-M0S2 can display 
even more interesting electronic and magnetic properties. 

We consider bare, as well as hydrogen saturated armchair and zigzag nanoribbons. These 
nanoribbons are specified by their width w in A or n number of M0-S2 basis in the unit cell. We 
take armchair nanoribbon with n=\2 and zigzag nanoribbons with n=6 as prototypes. The distance 
between Mo and S atom, cImo-s varies depending on the position in the ribbon. For example 
for n=\2, while at the center of the armchair nanoribbon, Jmo-5=2.42 A, and rfs_s=3.13 A, at 
the edge of the armchair nanoribbons, they change to Jmo-5=2.56 A and ds_s=3.27 A. The lattice 
parameters at the center of ribbons attain the same values as IH-M0S2. The average binding energy 
of hydrogen atoms passivating Mo and S atoms at the edges of the nanoribbon is £'^=3.64 eV. The 
lengths of Mo-H and S-H bonds are 1.70 A and 1.36 A, respectively. The distance between S 
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atoms at the edge is calculated as 3.27 A upon hydrogen passivation. 



Phonon calculations, stability and elastic properties 




Twisting mode 



DOS (arb. units) 



Figure 1: (Color online) Calculated phonon frequencies, Q(k) of the bare armchair M0S2 nanorib- 
bon with w= 17.75 A or n=\2 (there are 36 atoms in the primitive cell) are presented along symme- 
try directions of the Brillouin zone using Small Displacement Method (SDM), and corresponding 
densities of states (DOS). 



While the structure optimization through energy minimization yields an indication whether a 
given nanoribbon is stable, a rigorous test for the stability can be achieved through phonon cal- 
culations. If the calculated frequencies of specific phonon modes are imaginary, the structure 
becomes unstable for the corresponding k-wave vectors in the BZ. Here we present an example 
for the stability test of nanoribbons, whereby we calculate the phonon frequencies of the bare arm- 
chair nanoribbon with n=\2. The calculated phonon branches and corresponding density of states 



(DOS) are presented in [figure] [1][] 1 The out of plane (ZA) branch with parabolic dispersion and 



fourth acoustic branch (or twisting mode 61 ) characteristics of nanoribbons are obtained. Earlier, 
the branch of twisting mode was revealed in M0S2 nanotubes. 62 Similar twisting modes are also 
calculated for ZnO nanoribbons. 63 The overall shape of DOS of nanoribbons are similar to that 
of 2D M0S2 sheets, 64 except that the gap between optical and acoustical branches is reduced due 
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to edge phonon states. For the same reason additional peaks occur for fiat phonon branches of 
edge modes in band continua. All modes having positive frequency indicate that the bare armchair 
nanoribbon of M0S2 with n=12 is stable. It is also expected that other bare nanoribbons having 
n > 12 are stable. 

Having demonstrated the stability of nanoribbons, we next investigate their mechanical prop- 
erties by calculating the elastic properties. Currently, the behavior of honeycomb structures under 
tension has been a subject of current interest.— - — While IH-M0S2 have honeycomb structure, 
its atomic configuration and bonding of atoms are dramatically different from graphene. There- 
fore, the response of A- and Z-M0S2NRS to the strain is expected to be different from graphene. 
The elastic properties of the quasi ID M0S2 nanoribbons are examined through the variation of 
the total energy with applied strain. Generally, nanoribbons change their electronic and magnetic 
properties under uniaxial tension in the elastic deformation range. Here we present the response of 
A-M0S2NR and Z-M0S2NR to the strain in elastic range. 

Nanomechanics of both bare A-M0S2NR with n=\2 and Z-M0S2NR with n=6 is explored 
by calculating the mechanical properties as a response to the strain along the axis of the ribbon. 
To allow more variational freedom and reconstruction, segments of these NRs are treated within 
supercell geometry using periodic boundary conditions and spin-polarized calculations are carried 
out. Each supercell, both having total of 108 atoms, contains three unit cells for armchair and 
six unit cells for zigzag nanoribbons, respectively. The stretching of the ribbon is achieved by 
increasing the equilibrium lattice constant cq by Ac, to attain the axial strain e = Ac/cq. We 
optimized the atomic structure at each increment of the strain, Ae =0.01 and calculated the total 
energy under strain Et{c). Then the strain energy can be given by, E$ = Ej{e) —Ej{e = 0); 
namely, the total energy at a given strain e minus the total energy at zero strain. The tension force, 
Fj = —dEs(e)/dc and the force constant k = d 2 Es/dc 2 are obtained from the strain energy. 
Owing to ambiguities in defining the Young's modulus of honeycomb structures, one can use in- 
plane stiffness C = (I/Aq) ■ (d 2 Es/de 2 ) in terms of the equilibrium area of the supercell, Aq.— & 
The in-plane stiffness can be deduced from k by defining an effective width for the ribbon. 
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For both A- and Z-M0S2NR the hexagonal symmetry is disturbed, but overall honeycomb like 
structure is maintained in the elastic range. However, stretched ribbons can return to its original 
geometry when the tension is released. In the harmonic range the force constant is calculated to 
be K= 116.39 N/m and 92.38 N/m for A-MoS 2 having n =12 and Z-MoS 2 NR having n =6, re- 
spectively. Similarly, the calculated in-plane stiffness for the same ribbons are C=108.47 N/m and 
103.71 N/m, respectively. The difference between the values of armchair and zigzag nanoribbon 
occurs due to different bond and edge directions. As the width of the nanoribbon goes to infinity 
these two values are expected to converge to a single value. The calculated values are smaller than 
the values of C=292 and 239 N/m calculated for graphene and BN honeycomb armchair nanorib- 
bons. 70 Nevertheless, both calculated K and C values indicate the strength of IH-M0S2. It should 
be noted that K is approximately proportional to n, but C is independent of n for large n. Small 
deviations arise from the edge effects. 

For applied strains in the plastic deformation range the atomic structure of the ribbon under- 
goes irreversible structural changes, whereby uniform honeycomb structure is destroyed. At the 
first yielding point the strain energy drops suddenly, where the ribbons undergo an irreversible 
structural transformation. Beyond the yielding point the ribbons are recovered and started to de- 
form elastically until next yielding. Thus, variation of the total energy and atomic structure with 
stretching of nanoribbons exhibit sequential elastic and yielding stages. 

Electronic and magnetic properties 

Similar to the single layer IH-M0S2, its armchair nanoribbons (A-M0S2NR) are also semicon- 
ductors. The bare A-M0S2NR is a nonmagnetic, direct band gap semiconductor. Upon hydrogen 
termination of the edges, the band gap increases. Also the direct band gap shows variation with 
n, like the family behavior of graphene nanoribbons. However, unlike armchair graphene nanorib- 
bons, 73 the band gaps of A-MoS2NR's do not vary significantly with its width w or n. For narrow 
armchair nanoribbons with n < 7 the calculated value of the band gap is larger than that of wide 
nanoribbons due to quantum confinement effect. The variation of E g with n is in agreement with 
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that calculated by Li et.al.— 

The electronic band structure and charge density of specific states are examined in detail for a 
bare A-M0S2NR of n=12 in [figure][2][]2f a). The edge states, which are driven from Mo-4<i and 
S-3p orbitals and have their charge localized at the edges of the nanoribbon form flat bands located 
in the large band gap of IH-M0S2. Because of these edge states, the band gaps of bare armchair 
nanoribbons are smaller than that of IH-M0S2. Upon termination of each Mo atom at the edge by 
two hydrogen and each S atom by a single hydrogen atom, the part of edge states are discarded 



and thus the band gap slightly increases. As seen in [figure] [2] [] 2 (b), the remaining edge states 
continue to determine the band gap of the ribbon. Even if the character of these bands change, 
their charges continue to be located near the edge of the ribbon. Nevertheless, the band gaps of 
hydrogen saturated armchair nanoribbons remain to be smaller than that of 2D IH-M0S2. 

(b) 




E g =0.56eV 
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E g =0.61 eV 



Figure 2: (Color online) (a) Energy band structure of bare A-M0S2NR having n=12 and the width 
w=17.75 A. The band gap is shaded and the zero of energy is set at the Fermi level. At the right- 
hand side, charge density isosurfaces of specific states at the conduction and valence band edges 
are shown, (b) Same as (a), but the edge atoms are saturated by H atoms as described in the text. 
Large (purple), medium (yellow) and small (blue) balls are Mo, S, and H atoms, respectively. Short 
and dark arrows indicate the direction of the axes of nanoribbons. Total number of atoms in the 
unit cells are indicated. 



Furthermore we investigated the variation of band gap of hydrogen saturated armchair nanorib- 



bons as a function of n. As shown in [figure] [3] []3 for n < 7 the values of band gap are larger due 
to quantum confinement effect, but for n > 1 they tend to oscillate showing a family like behavior. 
These oscillations follow those found for bare armchair nanoribbon.— All calculated A-M0S2NR 
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are found to be direct band gap semiconductors. 



SO 0.50 




10 12 14 16 
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Figure 3: (Color online) Variation of direct band gap of hydrogen saturated armchair nanoribbons, 
A-M0S2NR with n. Each Mo atom at the edge is passivated by two hydrogen atoms and each S 
atom is passivated by a single H atom. 



In contrast to A-M0S2NR, the bare zigzag nanoribbons Z-M0S2NR are spin-polarized metals. 
Here we consider Z-M0S2NR with n=6 as a prototype. In [figure] [4][]4[ it is shown that the edge 
atoms of this nanoribbon undergo a (2 x 1) reconstruction by lowering its total energy by 0.75 
eV. Interestingly, as a result of reconstruction, the bare Z-M0S2NR is a half-metal with integer 
magnetic moment per primitive cell, namely ji=2 jiig. Thus the nanoribbon is metallic for ma- 
jority (spin-up) bands, but a semiconductor for minority (spin-down) bands with an indirect band 
gap of ~ 0.50 eV. We check that half-metallic state of bare Z-MoS2NR's is maintained for n=5, 
and n=S. Half-metals are interesting spintronic materials and were revealed first in 3D crystals. 74 
Lately, various nanostructures, such as Si nanowires 75 and atomic chains of carbon-transition met- 
als compounds 7 -- have found to display half-metallic properties. The half-metallic property is 
destroyed upon the saturation of the edge atoms by hydrogen. The magnetic moment of the ribbon 
and the density of spin states at the Fermi level depend on how Mo and S atoms at the edges of the 
ribbon are passivated by hydrogen. One distinguishes three different hydrogen passivations, each 



leads to different magnetic moments as indicated in [figure] [4] []4 As the number of passivating 



hydrogen atoms increases the number of bands crossing the Fermi level decreases. However the 
spin-polarization is relatively higher, if each S atoms at one edge are passivated by single hydrogen 
atom and each Mo atom at the other edge is passivated by double hydrogen. Interestingly, the lat- 



ter nanoribbon in [figure] [4] []4 d) is metallic for one spin direction and semimetal for the opposite 
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spin direction. Different spin polarizations found for different spin directions can make potential 
nanostructure for applications in spintronics. 

Z-MoS 2 NR 

(a) (b) (c) (d) 




(.1=2.00 |J=0.47 (1=0.82 (1=2.24 

Figure 4: (Color online) Atomic and energy band structure of bare and hydrogen saturated zigzag 
nanoribbon Z-M0S2NR having n=6 M0-S2 basis in the primitive unit cell. The top and side views 
of the atomic structure together with the difference of spin-up and spin-down charges, Ap = p^ — 
p^, are shown by yellow/light and turquoise/dark isosurfaces, respectively. The isosurface value 
is taken to be 10~ 3 electrons per A 3 . The (2x1) unit cell with the lattice constant 2a is delineated. 
Large (purple), medium (yellow) and small (blue) balls are Mo, S, and H atoms, respectively. The 
zero of energy is set at the Fermi level shown by dash-dotted green/dark lines. Energy bands with 
solid (blue) and dashed (red) lines show spin-up and spin-down states, respectively, (a) The bare 
Z-M0S2NR having ji=2 jig per cell displays half-metallic properties, (b) Spin-polarized ground 
state of Z-M0S2NR with Mo atoms at one edge and bottom S atoms at the other are passivated 
by single hydrogen, (c) Similar to (b), but Mo atoms are passivated by two hydrogen atoms, (d) 
Similar to (c), but top S atoms at the other edge are also passivated by single hydrogen atoms. The 
net magnetic moment of each case is indicated below the corresponding band panels. Bands are 
calculated using double cells. Small arrows along z-axis indicate the direction of the nanoribbon. 
The total number of atoms in supercell calculations are indicated for each case. 

Earlier Li et a/. 48 examined electronic and magnetic properties of armchair and zigzag M0S2 
nanoribbons using VASP 56 - 57 within DFT. They found that armchair nanoribbons are nonmagnetic 
semiconductors and their direct energy band gap vary with n and becomes 0.56 eV as n — > °°. They 
did not consider hydrogen passivation of edge atoms. They also noted that the value of net magnetic 
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moment can change, but the ferromagnetic state of zigzag nanoribbons are maintained even after H 
passivation of edge atoms. Mendez et ai— investigated armchair nanoribbons and concluded that 
these nanoribbons are metallic and have a net magnetic moment, but they change to semiconductor 
after hydrogen passivation of edge atoms. Their calculations show that in the case of bare armchair 
nanoribbons, the magnetic state is energetically more favorable by 14 meV and for H-saturated 
zigzag nanoribbons the antiferromagnetic state is favorable relative to the ferromagnetic state by 
15 meV. These results disagree the present results, as well as with those of Li et ah— 

Normally, the bare and unreconstructed zigzag nanoribbons have sizable electric dipole mo- 
ment along the direction from the edge having only negatively charged S atoms to the other edge 
having only positively charged Mo atoms. The dipole moment is calculated to be 55.4 eA per cell 
of Z-M0S2NR having n = 6, but it reduces to 0.07 A upon reconstruction of the edges. Present 
results show that the edge reconstruction ought to be treated properly to reveal the half-metallic 
state and to estimate the correct dipole moment. 

Impurities and Defects in M0S2 Nanoribbons 

Interesting properties of M0S2 nanoribbons revealed above can be modified through adatom ad- 
sorption (or doping) and vacancy defect creation. Earlier, Huang and Cho 4 -^ investigated the ad- 
sorption of CO on a pure IH-M0S2 surface by using DFT Similarly, aromatic and conjugated 
compounds on M0S2 are also studied. 43 Similar to graphene-^ - — and its nanoribbons,— - — not 
only adatoms, but also vacancy defects created can led to crucial effects. 

Here we consider again our prototype armchair nanoribbon and investigate the adsorption of 
C, O and Co. C is widely investigated in other honeycomb structures; the adsorption O is ex- 
pected to result in important changes due to its high electronegativity. On the other hand, Co 
being a transition metal atom is expected to attribute magnetic properties. In the supercell ge- 
ometry, a single adatom is adsorbed at every three unit cells, which leads to the adatom-adatom 
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Figure 5: (Color online) Top and side views for the schematic representation of possible adsorption 
geometries of adatoms obtained after the structure optimization. Adatoms, Mo and S are repre- 
sented by red (large-dark), purple (medium-gray) and yellow (small-light) balls, respectively. Side 
view clarifies the height of adatoms from Mo and S atomic planes. In each possible adsorption 
geometry, the entry on the lower-left part indicates where the adatom is initially placed. All sites 
show geometries associated with the adsorption to a bare armchair (n=12) nanoribbon (NR). The 
calculations are carried out in the supercell geometry where a single adatom is adsorbed at every 
three unit cells. The total number of atoms in the supercell is 109. Possible adsorption geometries 
in NR/£ (adatom is initially placed at the inner edge of bare armchair NR) and NR^ (adatom 
is initially placed at the outer edge of bare armchair NR). Adatoms indicated at lower right part 
of every possible adsorption geometry correspond to those, which are relaxed to this particular 
geometry upon structure optimization. 
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Table 1: Calculated values of adatoms adsorbed to the bare armchair M0S2 nanoribbon 
having n=12 M0S2 units in the primitive unit cell. The supercell in calculations consist of 
three primitive cells. There are two different adsorption sites as (described in |[figure][5][]5| ) 



for each adatoms. The positions only with a positive binding energy is indicated. h Mo -, the 
height of the adatom from Mo layer; %, the height of the adatom from the nearest S-layer; 
dMo-, the adatom-nearest Mo distance; ds, the adatom-nearest S distance; E b , adatom binding 
energy; /i T , magnetic moment per supercell in Bohr magneton ji B ; p% excess charge on the 
adatom (where negative sign indicates excess electrons); 4>, photoelectric threshold (work 
function); P, dipole moment calculated in the x, y and z direction, respectively. Nanoribbon 
is in the (x,y) -plane and along the x-direction. E t , energies of localized states induced by 
adatoms. Localized states are measured from the top of the valence bands in electron volt. 
The occupied ones are indicated by bold numerals and their spin alignments are denoted by 
either | or \.. States without indicated spin alignment are nonmagnetic. 
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distance of ~ 16.60 A. We found that the edges of the nanoribbon are active sites for adsorp- 
tion and are energetically more favorable relative to the center of nanoribbon. As described in 



[figure] [5] [] 5 adatoms adsorbed at the inner (NR/^) and outer parts of the edges (NR^e) of the 



armchair nanoribbon result in a reconstruction on the edges and form strong bonds with nanorib- 



bon. In | [table] [ 1 ] [JT] we present all relevant data obtained from our calculations of adatoms ad- 
sorbed to A-M0S2NR. The height of the adatom from the Mo- or S-planes are calculated relative 
to the average heights of Mo- and S- atoms in the corresponding planes. The binding energy, 
Ef, is calculated as Eb = E a d + Ea-MoS 2 ~ Ead+A-MoS 2 - Here, E ac { is the ground state energy of 
free adatom calculated in the same supercell with the same parameters; Ea~moS 2 l% me total en- 
ergy of nanoribbon and E ac / + A-MoS 2 is me energy of adatom+A-MoS2NR complex. The charge at 
the adatom pg, is calculated using Bader analysis. 58 The excess charge of the adatom is obtained 
by subtracting the charge at the adatom, pg from the valence charge of the adatom Za, namely 
p* = Za~Pb- Accordingly p* < implies excess electron at the adatom site. Here the adatom+A- 
M0S2NR complex attains net magnetic moment after the adsorption of transition metal atom, Co. 
Adsorption of C and O do not cause any spin polarization in all adsorption geometries. Adsorbed O 
having the highest electronegativity among the adsorbates treated here has highest excess charge; 
C is also negatively charged in both adsorption geometries. Co adatom having electronegativity 
smaller than those of both constituent atoms of the nanoribbon is positively charged. The deple- 
tion and annihilation of charge from the adatom result in a small dipole moment in the y-direction, 
which is normal to the ribbon. Since adatom-adatom interaction is hindered due to large supercell 
dimensions, the localized states form flat bands in the supercell geometry. For C and O localized 
states deep in the valance band are generally due to their low energy 2s-states. For Co adatom most 
of the localized states originate from 3J-orbitals. 

We note that previously, He— et al. found that the lowest energy adsorption position of C 
adatom is at the top of Mo atom in monolayer M0S2 and oxygen adatom is adsorbed to the top site 
above the S atom. However, adsorption of adatoms at the edges of the M0S2 nanoribbon gives rise 
to properties rather different from those in IH-M0S2. 
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Vacancy Defects in M0S2NR 



We investigated five different types of vacancy defects, namely Mo-, S-vacancy, MoS-, S2-divacancy 
and MoS2-triple vacancy, which are formed near the center of hydrogen passivated armchair 
(n = 12) and zigzag (n = 6) nanoribbons. All structures are optimized after the creation of each 
type of vacancy. Vacancy energies Ey, are calculated by subtracting the sum of the total energy 
of a structure having a particular vacancy type and the total energy(ies) of missing atoms in the 
vacancy defect from the total energy of the perfect structure (without vacancy). Here all structures 
are optimized in their ground states (whether magnetic or nonmagnetic). Positive Ey indicates that 



the formation of vacancy defect is an endothermic process. In |[table][2][]2| calculated vacancy 
energies as defined above and their magnetic ground states are presented. 

Table 2: Calculated vacancy energies Ey (in eV), magnetic moments /1 (in ji B ) of five different 
types of vacancy defects, Mo, MoS, MoS 2 , S, S 2 in A-M0S2NR and Z-M0S2NR. NM stands 
for nonmagnetic state with net /i=0. E h energies of localized states in the band gap. Local- 
ized states are measured from the top of the valence bands in electron volt. The occupied 
ones are indicated by bold numerals and their spin alignments are denoted by either t or |. 
Nonmagnetic states have no spin alignments. 



Mo MoS MoS 2 S S 2 

E v -\i E v -ll Ey-\i E\/-ll Ey-[L 

A-MoSiNR 16.92-NM 17.47-NM 22.94-2.00 5.82-NM 11.55-NM 

Ej 0.09. 0.11. 0.35. 0.48 0.11, 0.40, 0.49 0.02t, 0.034. 0-33;, 0.34|, 0.50t - - 

Z-MoS 2 NR 15.78-8.06 16.41-8.66 22.02-8.67 5.09-8.61 10.77-8.31 



For a hydrogen saturated armchair nanoribbon (A-MoS 2 NR), having width n=\2, the vacancy 
defects are treated in a supercell geometry, where a single defect is repeated in every four unit 
cell. For this supercell configuration, vacancy-vacancy coupling becomes minute and the resulting 
defect states appear as flat bands. Similar to 1H-MoS 2 , all the vacancy types have zero net magnetic 
moments, except MoS 2 -triple vacancy, which has a net magnetic moment of /! = 2/ig per supercell. 
The nonmagnetic excited states associated with vacancy defects occur above ~120 meV, and are 
derived from Mo-4J and S-3p orbitals around the vacancy. 

Similar to the armchair nanoribbons, various vacancy and divacancy defects in the hydro- 
genated zigzag nanoribbon (Z-MoS 2 NR), having width n=6 are treated in a supercell geometry, 

15 



where a single defect is repeated in every eight unit cell. Calculated vacancy energies, net mag- 



netic moments per cell are presented in |[table][2][]2| It is found that in the presence of a vacancy 
defect, such as MoS-divacancy, S2-divacancy, Mo-vacancy and S-vacancy, the spin-polarization of 
the zigzag nanoribbons appears to be suppressed. For example, while defect free, zigzag nanorib- 
bons are metallic and spin-polarized ground state with net magnetic moment of /i = 2.24/ig per 
double unit cell, the total magnetic moment of eight unit cell decreases to ~ /i = 8.30/lg per super- 
cell from fx = 8.96/ig in the presence of vacancy defects. In particular, the net magnetic moment 
of MoS2-triple vacancy appear to compensate for the edge magnetization of the zigzag M0S2 NR 
to result in a net magnetic moment of 8.67 /ig per supercell. 



Discussion and Conclusions 

The phonon dispersion of bare A-M0S2NR with n = 12 is calculated and stability of nanoribbons 
are ensured. Armchair nanoribbons are nonmagnetic direct band gap semiconductors; their en- 
ergy band gaps vary with its width and termination of edge atoms with hydrogen, whereas zigzag 
nanoribbons are found to be ferromagnetic metals. The bare zigzag nanoribbon is found to be a 
half-metal. Both nanoribbons are stiff materials, but their in-plane stiffness are calculated to be 
less than half of those of graphene and BN. 

The adsorption of adatoms and creation of vacancy defects in M0S2 nanoribbons have crucial 
effects in the electronic and magnetic properties. We found that several adatoms can be adsorbed 
readily at diverse sites with significant binding energy. In this respect, M0S2 appears to be a 
material, which is suitable for functionalization. Similarly, net magnetic moment can be achieved 
through the adsorption of Co adatoms to the nonmagnetic armchair nanoribbons. In addition to 
spin-polarization, significant charges are transferred to (or from) adatom. 

While vacancy defects of S, S2, Mo and MoS created in hydrogen passivated armchair nanorib- 
bon do not induce any magnetic moment, the creation of M0S2 triple vacancy results in a signifi- 
cant magnetic moment in the system. Vacancy creation in hydrogen passivated zigzag nanoribbons 
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however suppresses the magnetic moment occurring at the edges of the nanoribbon and results in a 
decrease in the total magnetic moment of the system. Briefly, functionalization of M0S2 nanorib- 
bons through adatom adsorption and vacancy creation appears to be a promising way to extend the 
application of M0S2. 
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